Quantum noise and mixedness of a pumped dissipative non-linear oscillator 
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Evolutions of quantum noise, characterized by quadrature squeezing parameter and Fano factor, 
and of mixedness, quantified by quantum von Neumann and linear entropies, of a pumped dissipative 
non-linear oscillator are studied. The model can describe a signal mode interacting with a thermal 
reservoir in a parametrically pumped cavity with a Kerr non-linearity. It is discussed that the initial 
pure states, including coherent states, Fock states, and finite superpositions of coherent states evolve 
into the same steady mixed state as verified by the quantum relative entropy and the Bures metric. 
It is shown analytically and verified numerically that the steady state can be well approximated by 
a nonclassical Gaussian state exhibiting quadrature squeezing and sub-Poissonian statistics for the 
cold thermal reservoir. A rapid increase is found in the mixedness, especially for the initial Fock 
states and superpositions of coherent states, during a very short time interval, and then for longer 
evolution times a decrease in the mixedness to the same, for all the initial states, and relatively low 
value of the nonclassical Gaussian state. 

Keywords: quantum entropy, quadrature squeezing, sub-Poissonian statistics, Wigner function, 
Husimi functions, Kerr non-linearity, Schrodinger cats, steady states. 



I. INTRODUCTION 



Quantum noise and mixedness are properties central 
for quantum theory |l|. It is obvious that a pure signal 
state interacting with a thermal environment loses its 
purity and turns into a mixed state. Since real optical 
devices suffer from losses, it is important to study their 
influence on the noise and mixedness during the state 
evolution. 

In the case of interactions of a small number of the 
signal and reservoir modes, the noise and mixedness pa- 
rameters periodically rise and fall. But in the case of the 
signal interacting with an infinite number of the reservoir 
modes, the evolution is irreversible and after a few char- 
acteristic time intervals, the signal mode transforms into 
an asymptotic steady state depending on external param- 
eters of the system, e.g., an external classical pump. One 
can show that when the pump is turned off, final steady 
state of the signal will be the pure vacuum state, but 
when the pump is turned on, the steady state is not the 
vacuum but a mixed nonclassical state. We show further 
that the steady state can be properly approximated by a 
nonclassical Gaussian state (NCGS). 

In the next sections we will investigate a particular 
example of a non-linear interaction, i.e., a nontrivial dy- 
namics of initial signal modes in a pumped resonator with 
a non-linear Kerr medium when losses are included (see, 
e -g-> H, H, 3). We are interested mainly in the quan- 
tum noise and mixedness (purity) evolutions of the sig- 
nal mode. We discuss three different types of the signal 
initial states including coherent states, Fock states and 
superpositions of coherent states. We study the influ- 
ence of the losses, pump and non-linearity on the signal 
evolution, in particular the fast decoherence of the initial 
coherent superpositions. 



II. INTERACTION MODEL 

A pumped non-linear oscillator can be described in the 
interaction picture by the following Hamiltonian 0, : 



H — Hj + Hr, 
Hi = i(pa t -p*a) + Ga t2 a 2 . 



(1) 
(2) 



where fi and fit arc the annihilation and creation oper- 
ators of the signal mode, respectively. The first term 
in Hamiltonian Hi describes a pump process with the 
complex amplitude p, the second term represents a Kerr 
process with the non-linear interaction constant G and 
Hr is the term describing the interaction of the signal 
mode with the cold thermal reservoir. For simplicity we 
use the units where h = 1. The model can be used in a 
description of, e.g., a parametrically pumped cavity with 
a Kerr non-linearity. We adopt the Heisenberg-Langevin 
approach in which the Heisenberg equation for the 
signal mode reads is 



— = p — 2iGfi^fi 2 — 7nfi 
dt 



C. 



(3) 



where C is the Langevin force and 70 is the loss pa- 
rameter. For the optical modes at room temperature, 
the number of thermal photons can be assumed negli- 
gible. Therefore, in the following we assume the zero- 
temperature approximation. In that case the standard 
rules hold for the Langevin force: 



(Aii)A*a)> = (CHh)C(t 2 )) = (£(*!)) = 0, 
(C(h)C\t 2 )) = 2 7o S(t 1 -t 2 ). 



(4) 



For the linear interaction (G = 0) of the signal mode, 
an analytical solution can be found. For the coherent 
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initial state the solution is 
time-dependent complex amplitude 



\a(t)) with the 



a{t) 



-lot 



^(1 

7o 



-lot\ 



(5) 



So the initial coherent state with the amplitude a finally 
evolves into another coherent state with the amplitude 
p/jo- The evolution is essentially finished during a short 
period equal to the first few characteristic times r = 
I/70. For the model without a pump, the initial state 
turns into the vacuum state. 

In the non-linear Kerr interaction case, the evolu- 
tion cannot be solved analytically for arbitrary evolution 
times and suitable numerical methods have to be used. 
The state evolution in this case is much more complex 
and is the main subject of our study. Under reason- 
able assumptions some analytical approximations can be 
derived as well. After presenting our precise numerical 
results, we will also give an exact quantum steady-state 
solution and its semiclassical approximation in the long- 
time limit. 

Alternatively we can start from the Liouville equation, 
which leads to the following master equation 



(6) 



where R for the zero temperature reservoir stands for 

R = — [Hr,p] — 70 (2apa) — a)ap — pa^a) . (7) 

Master equation © can simply be expressed in the Fock 
basis with the density matrix in the form 



9 



rim 



\m)(n\ 



(8) 



leading to a set of ordinary linear differential equations 
for the density matrix elements p mn = (m\p\n). For weak 
interaction fields the differential equations can be solved 
numerically. The exact quantum results presented in this 
article are obtained by applying this method. 

To visualize evolution of quantum states (pure or 
mixed) generated in our system, we apply the Husimi 
and Wigner functions, which are the special cases of the 
Cahill-Glauber s-parametrized quasidistribution func- 
tion W [s) {f3) defined for -1 < s < 1 as follows @: 



1 



F'M = -Tr{pTW(f3)}, 

IT 

where p is the density matrix of the field and 
1 



(9) 



fW(/3) 



exp(/?r-/n)£ (s) (£)d 2 £, (10) 



and £>M(£) = c s ^ 2 / 2 D(£) with £>(£) being the displace- 
ment operator. The quasidistribution M/"()3) in the 
number-state basis can be calculated as 0: 



W (s) (/?) = i$> m „HfM(/?)|m>, 



(11) 



where 

{n\T^{f3)\ m )={— { 



1 \ V 2 / 2 \ m - n+1 / a + 1 



1 - s 



, i r )"'-", -x]» ( -^3^ 1 1"." " 



8-1 

4|/?| 2 
1-s 5 



(12) 



given in terms of the associate Laguerre polynomials 
U™~ n {x). Eq. (|T2*|) for s — > — 1 goes into the simple 
expression 



(n\T 



(-1) 



exp l 



\J n\m\ 



(13) 



as can be derived by observing that lim e ^o e ™i™" n (7) = 
{—y) n /n\. Eq. (|T2l can also be applied for s — ► 1 
if the limit is taken carefully (see e.g. 0). The spe- 
cial cases of W (s) (/3) for s = -1,0,1 arc known as the 
Husimi Q-function, the Wigner VF-function and Glauber- 
Sudarshan P-function, respectively. For example, the 
s-parametrized quasidistribution function for coherent 
state I a) is given by the Gaussian distribution yV^ s )(/3) = 
2/[vr(l - s)] exp[-2|^ - a| 2 /(l - s)], which for s = 1 be- 
comes Dirac's delta 8{(5 — a). If a given state is described 
by the P-function, which is positive definite and no more 
irregular than Dirac's delta, then the state is classical, 
otherwise the state is considered to be nonclassical. We 
use this criterion to show that the steady state of our 
system is nonclassical. 



III. MEASURES OF QUANTUM NOISE AND 
MIXEDNESS 



Quantum noise properties of nonclassical light can be 
analysed in terms of the Fano factor and quadrature 
noise variances. The first parameter corresponds to di- 
rect photo-pulse detections, while the second is related to 
homodyne detection schemes. The Fano factor is defined 
by 



F = 



(n) 



(14) 



which for Poissonian states satisfies F = 1, while for sub- 
and super-Poissonian states is F < 1 and F > 1, respec- 
tively. In particular, the Fock states are sub-Poissonian 
with F = independently of the number of photons, 
coherent states are Poissonian (F = 1), while thermal 
chaotic states are super-Poissonian since F = 1 + 2(n) > 
1. The quadrature noise squeezing parameter S is de- 
fined as the minimum variance S = ming((AAg) 2 ) over 
all possible values of phase 9 e (0, 2n) of the general 
quadrature operator Xg — ae~ lB + a^e l6 . For a coherent 
state 5=1, while a state with S < 1 is referred to as the 
quadrature squeezed light, since it has lower noise level 
than the coherent or vacuum state. 
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The most natural measure of mixedness of a state, 
given by p, is the von Neumann entropy 



E = -(lap) = -Tr{pln/S}. 



(15) 



Density matrix p of any mixed state can be expressed as 
the incoherent sum 



P 



N 

£ 

k=l 



Pk\tpk)(lpk 



(16) 



of the orthogonal pure states \ipk), where pk are their 
weight factors being eigenvalues of the density matrix 
p. As follows from the general properties of the density 



matrix, it holds < pk < 1 and X^feLiPfc = 1- Thus, the 
von Neumann entropy can be expressed as 

JV 

E = -J^Pklnpk > 0. 



(17) 



k=l 



Another useful measure of mixedness is the linear entropy 
L defined as 



L = 1 - Tr{p 2 }, 



(18) 



where the second term in l|18[) is referred to as the purity 
P of the state 



JV 



P = Tr{p 2 } = £ 



(19) 



fe=i 



The mixedness L and the purity P are complementary 
in the sense that whenever the entropy increases the 
purity falls. For completeness, we note other general- 
ized measures of the mixedness including (see, e.g., ||): 
the Renyi entropies defined by H q — (lnTrjp 9 })/ (1 — q) 
{q = 2,3,...), the Neumann-Renyi entropy given by 
H-2 = — lnTr{/5 2 } = — InP or the participation ratio 
defined to be R = 1/Tr{p 2 } = l/P. 

It is easy to show the following properties of the mixed- 
ness parameters E and L: For any pure state it holds 
E = L = 0. For a mixture of two orthogonal states 
it holds < L < 1/2 and < E < ln2. For a bal- 
anced mixture of two orthogonal modes pi = p% = 1/2 
it holds L = 1/2 and E = In 2 as their maxima. For 
a mixed state composed of N orthogonal pure states it 
holds < L < (N - 1)/N and < E < \nN. For a 
homogeneous superposition it is pk = 1/N and therefore 
L = (N — 1)/N and E = In TV. So, for a state strongly 
mixed with the reservoir, L can increase to one and E to 
infinity. 

For thermal or chaotic states, the density matrix is 
diagonal with 



(20) 



Pkk =Pk = — 

(l + (n)) 

implying that the linear entropy is 

2(n) 



i+fc : 



^cha 



l + 2(n 



(21) 



and the corresponding von Neumann entropy is 



-Echaot — — In 



)<"> 



(l + (n» 1+<n> ' 



(22) 



Here, (n) stands for the number of photons in the chaotic 
mode. For the vacuum state (n) = 0, which is a pure 
state, we get E = L = 0, while for the chaotic state with 
(n) = 1 we get E = ln4 rs 1. 386 and L = 2/3. 

The thermal state is defined as the state, which max- 
imizes the entropy when at the same time the energy is 
fixed. So, the state with a fixed number of photons (n) 
has the upper limit of entropy and the following inequal- 
ity holds E < -E max = -Echaot • The linear entropy has the 
upper limit of 



1 



l + 2(n) 



(l + 3(n))(l + 3(n)/2) 



(23) 



reached for a state with the descending arithmetic se- 
quence type distribution 



Pk 



2 + 3 (n) 



1 - 



1 + 3(ra 



(24) 



For example, the state with (n) 
of the linear entropy equal to L B 



1 has the upper limit 
< = 7/10. 



IV. NUMERICAL ANALYSIS OF MIXEDNESS 
AND NOISE 

A. Ideal non-linear oscillator without pump 

It is well known that in the case of Kerr dynamics 
without losses and without a pump, the initial coherent 
state evolves periodically into a nonclassical light with 
highly reduced quantum noise (see also |10| and ref- 
erences therein) and also becomes superpositions of two 
(Schrodinger cats) [ll| and more (Schrodinger kittens) 
[l2| coherent states. If the initial state is 



k=0 



then its evolution is described by 



c k e 



-ik(k-l)Gt 



\k). 



(25) 



(26) 



fc=0 



being clearly periodic with the time period T = tt/G, as 
the state is \tp(T)) — \ipo). For example, the coherent 
initial state \ipo) = \a) at the time t — T/2 becomes the 
coherent superposition 



e i7r/4 



ia) + c 



— i7T / 4 I 



(27) 



More generally, at the evolution times t = —T, which 
is a rational fraction of the period T, the output state 
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FIG. 1: Time snapshots of Husimi Q-function for the initial 
coherent state \a) = |3). Kerr parameter G = 1, no loss 
70 = and no pump p = 0. 




FIG. 2: Same as in figure 1 but for the dissipative model with 
7o = 1. 



is given as a coherent superposition of n coherent states 
dislocated regularly on the circumference of a circle with 
a radius |a|. A typical example of such a Kerr evolution 
can be seen in figure 1. Obviously, the initial pure state 
remains pure (E = L = 0). 



B. Dissipative non-linear oscillator without pump 

When losses are involved in the system without a 
pump, the evolution is no more periodic and it ends in 
the vacuum state. During the time evolution, an ini- 
tial state goes through various nonclassical states and 
through strongly mixed states. Some analytical solutions 
of the dissipative Kerr non- linear oscillator have been ob- 
tained both for 'quiet' (T w 0) and 'noisy' (T > 0) 
[l4| reservoirs. 

The role of dissipation is clearly seen by comparing 
the evolution of the initial coherent state \a = 3) without 
(figure 1) and with (figure 2) losses. The parameters used 
in figures 1 and 2 are G = 1 , p = and the loss parameter 
7o is in figure 1 and 0.1 in figure 2. Due to losses, no 
superpositions of coherent states arise in the evolution 
presented in figure 2. The effect can be explained as a fast 
destruction of internal coherence. Now we will analyse 
the evolution of the mixedness, characterized by the von 
Neumann (E) and linear (L) entropies, and the noise 
in terms of the Fano factor F and quadrature squeezing 
parameter S given in figure 3 for the parameters G = 0.2, 
p = and 7o = 1. First we focus on the evolution of the 
initial coherent state |a) = |3) depicted by thick solid 
curves in figure 3. As seen, the signal state is maximally 
mixed at the time t « 0.87 when the entropies reach the 
values of E w 0.57 [figure 3(a)] and L w 0.31 [figure 
3(b)]. During the evolution, the photocount statistics 
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FIG. 3: Evolution of the mixedness and noise for the system 
without a pump: (a) the von Neumann entropy E, (b) linear 
entropy L, (c) Fano factor F and (d) squeezing parameter 
S for the initial coherent state \a — 3} (thick solid curves), 
the Fock state \n — 9} (dashed curves) and the superposition 
of three coherent states \1jj3) (thin solid curves). Resonator 
parameters are p = 0, G = 0.2 and 70 = 1. 




FIG. 4: Evolution of the purity P and entropy E for the initial 
Fock state \n = 9) (dotted curve) and for the superposition of 
three coherent states \tj)3.) with Kerr (thick solid) and without 
Kerr (thin solid curves). The resonator parameters are same 
as in figure 3. 



of the initial coherent state remains Poissonian. It is a 
direct consequence of the fact that the diagonal terms of 
the density matrix keep the Poissonian character for all 
the times. It can be shown that it holds exactly that 



p(n) = exp(-(ra» 



(28) 



where 



|a| 2 exp (— 2 jot), which explains why the 



Fano factor remains constant and equal to one for all the 
times [see figure 3(c)]. On the other hand, the quadra- 
ture noise evolves leading to the maximum squeezing of 
S w 0.52 at the time t w 0.18. In figure 3, we have 
also presented by thin solid curves the evolutions of the 
mixedness and the noise of the initial superposition of 
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three coherent states (Schrodinger kittens) 

m = N 3 (|3) + |3e^) + |3e-^>) , (29) 

where A3 is a normalization constant. One can observe 
a rapid increase in the mixedness measured by both E 
in figure 3(a) and L in figure 3(b) during a very short 
time interval t < 0.05. It can be explained as a general 
effect of the coherence loss of quantum components of the 
superposition In general, it can be shown that the 

superposition of three strong coherent states (normalized 
by A« 1/V3): 

\1>)=NQa l ) + \a 2 ) + \a 3 )) (30) 

exponentially fast evolves into the mixed state 

p = N 2 (K)(ai| + |a a ><o 2 | + |aa)<Oa|) (31) 

as a consequence of the reservoir influence. The charac- 
teristic time of the decoherence process between \cti} and 
\&k) (i ^ k f° r i,k = 1, 2, 3) can be estimated as 

T ik « (27o|«i - afc| 2 ) . (32) 

So for strong fields, the decoherence process is much 
faster than the dissipative process itself. The mixedness 
measured by the linear entropy increases to a value of 
L = 2/3, while by the von Neumann entropy increases to 
E = In 3 1.10. In our case of ^3) the increase in the 
linear entropy can be approximated by 

L« I [l-exp(-6 7o |a| 2 t)] . (33) 

whose validity is verified by the results shown in figures 
3(a,b) and 4. We note that a better agreement with 
the decoherence theory is achieved for the system when 
the Kerr non-linearity is switched off (corresponding to 
thin solid curves in figure 4) rather than is switched on 
(depicted by thick solid curves). For longer times, the 
noise dominates the process and the signal state turns 
into the vacuum state. Thus, the squeezing parameter 
evolves into the value of one as shown in figure 3(d). 
As for the initial coherent state, the Fano factor for the 
initial superposition state remains unchanged and 
equal to unity [see figure 3(c)] during the evolution of the 
unpumped system for arbitrary values of the interaction 
G and loss 70 parameters. 

For the initial Fock states \n = 9), a rapid increase 
in the mixedness can also be observed in figure 3(a,b) 
as depicted by dotted curves. For t ~ 0.35 we find the 
mixedness parameters to be L « 0.8 and E w 1.8. The 
mechanism of decoherence and increase in the mixedness 
differ from those for the initial superposition of coherent 
states. We can explain the observed values simply as 
follows: by solving directly the Schrodinger equation for 
the initial Fock state |n), we get a diagonal density matrix 
at time t described by the binomial distribution 

P(k) = Pkk = (f) (c-^'f (1 - e- 2 ^p fc , (34) 



which is surprisingly completely independent of the Kerr 
non-linearity G. The mixedness parameter L = 1 — 
^2p 2 (k) sums up to an expression proportional to a hy- 
pcrgeomctric function 2 -Pi {—n,—n, 1, x), which is diffi- 
cult to handle. But we can easily calculate the variance 

A2fc = n (e" 270 ') (1 - e- 270t ) (35) 

which reaches a maximum at the time moment when 
e -27ot _ \J2. So we can estimate the searched maxi- 
mum of the linear entropy L at this moment. We have 

Efn\ 2 <>„ (2n)! 
fe=0 U 2 =1 -(kiF w 

and for large n we get the following simple approximation 

imax ~ 1 (37) 

by applying the Stirling formula nl w \f2wne~ n n n . Thus, 
we can analytically estimate the loss of purity for n — 9 as 
L max w 0.815 according to l(3*B|) (or slightly less accurate 
as 0.812 according to J57JI), which can be expected at 
the time t w 0.347. We can calculate the entropy of 
the state at the same time moment to get E w 1.823. 
These estimations are in full agreement with the precise 
numerical results given in figure 5(a,b). Finally, we note 
that the squeezing parameter decreases from S = 19, 
while the Fano factor increases from F = for the initial 
Fock state to reach the value of one for the vacuum state 
in the time limit. 



C. Dissipative non-linear oscillator with pump 

Here, we will analyse the system with the classical 
pump switched on, for which the initial state finally 
evolves into a stationary non-vacuum state. Its proper- 
ties are determined by the pump intensity p and by the 
passive parameters 70 and G. We will study three typi- 
cal evolutions for different initial fields including coherent 
states, Fock states and superpositions of three coherent 
states. Figure 5 shows the Wigner function evolution of 
the coherent initial state \a — 3), figure 6 shows the evo- 
lution of the Fock initial state \n = 9) and figure 7 shows 
evolution of the superposition of three coherent states, 
given by (|29|l . for the same resonator parameters p = 5, 
70 = 1 and G — 0.2. All the three states have been se- 
lected to have the same initial energy (n) — 9 (in the last 
case approximately). Six time snapshots of the Wigner 
function of the signal state are presented for the time 
moments t — 0,1/3,2/3,1,5 and for a very long time, 
practically corresponding to 00. In figure 5, the coherent 
initial state rotates and evolves in a banana-shape state 
and finally ends in a steady state. In figure 6, the ini- 
tial Fock state of a ring shape is deformed and squeezes 
before ending in a steady state. In figure 7, the initial 
superposition of the coherent states rotates, deforms and 
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FIG. 5: Time snapshots of Wigner function for the initial 
coherent state \a = 3). Resonator parameters are G = 0.2, 
70 = 1 and p — 5. 



FIG. 7: Same as in figure 5 but for the initial coherent state 
superposition l^}- 





FIG. 6: Same as in figure 5 but for the initial Fock state 
|n = 9>. 



FIG. 8: Evolution of the mixedness and noise parameters 
same as in figure 4 but for the pumped system with p = 5. 



its components converge into the same steady state. 
For the coherent initial state, the phase-space evolutions 
of the mean quantum (a(t)) and classical a(t) amplitudes 
in the time interval < t < 10 are presented in figure 
9. Both the spirals start at the point a = 3, and end at 
a w 1 — 2i. Note that the evolutions are very similar to 
each other even for the other times, which shows the va- 
lidity of our semiclassical approximation. The evolutions 
of the mixedness parameters E and L and the quantum 
noise parameters F and S for all the three initial states in 
the pumped system are depicted in figure 8. Note a rapid 
increase in the mixedness especially in the cases of the 
initial Fock state and a superposition of coherent states 
for very short evolution times, and then for longer times, 
a decrease in the mixedness for all the initial states. We 
observe that the evolutions for short times are similar 



to those presented in figure 3 for the unpumped system. 
However, for longer times (even for t > 5 for the parame- 
ters of figure 8) all the three initial states of the pumped 
system evolve into some mixed steady state different from 
the vacuum state. This is in contrast to the evolution of 
the unpumped system for which E(p = 0) and L(p = 0) 
approach zero, while S(p — 0) and F(p = 0) approach 
one in the time limit (see figure 3), since the dissipa- 
tive system without external pumping evolves into the 
vacuum state for any initial fields. It is worth noting, 
although it is out of the focus of our interest, that the 
weakly-pumped non-linear oscillator (j> <C G) enables the 
so-called optical state truncation [15( and can be used for 
optical qubit generation. 
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Wi 




FIG. 11: Wigner functions W 3 = W s ( s 0) (|^)) (J = 0,1,2) of 
the first three pure state components of the steady state for 
the resonator parameters G = 0.2, p = 5, 70 = 1. 



Re a 



FIG. 9: Phase space trajectory of quantum (a (t)) (thick 
curve) and classical a (t) (thin curve) solutions for t £ (0, 10). 




FIG. 10: Evolution of the Bures metric (Db) and the quan- 
tum relative entropy (Dkl) for the same initial states as in 
figure 8. 



D. Steady-state solution 

In figures 5-7 we have presented the evolutions of the 
Wigner functions up to the time t = 5 and then checked 
for much longer times, practically, for t — > 00. From this 
study we can conclude that all the three initial states end 
in the same steady state, which is centred around the 
same value of the complex amplitude a w 1 — 2i. All the 
other parameters studied that is: E w 0.278, L w 0.135, 
(ri) « 5.13, F ra 0.69 and S ~ 0.72 of the steady state are 
the same as well. So, the steady-state solution is appar- 
ently unique and totally independent of the initial states. 
It can be therefore fully determined by the resonator pa- 
rameters only. While for the system without a pump the 
evolution ends in the vacuum state, for the pumped sys- 
tem the asymptotic state is neither the vacuum nor pure, 
and can have intriguing noise properties. We can prove 
it by a decomposition of the steady-state density matrix 
in its spectrum of eigenstates. 

Drummond and Walls found that no Glauber- 
Sudarshan P-function exists in the steady-state except 
as a generalized function 0. The latter was used by 
Khcruntsyan to find the following explicit form of the 
steady-state density matrix Q 



(e*r £ " F 2 (A*+m,A + n,|e| 2 ) 
V^W r(A* +m)T(X + n) 



(38) 

where C = T(\*)Y{\) / a F 2 {\* , A, 2|e| 2 ), e = -ip/G, 
A = — «7o/G, and qF 2 is the generalized Gauss hyper- 
geometric function. Eq. I|38|l readily enables calculation 



of the moments {(a' 
up to 



by the expressions summing 



, f -Um~n\ rU *\m n oF 2 (A* + 777, A + 77, 2|e| 2 ) 

((a } a >ss = c{£) £ r(A* + m)r(A + 77) (39) 

which corresponds to a slightly modified Drummond- 
Gardiner formula Solution with the help of 

definition (|ll|l . also enables calculation of the steady- 
state Husimi and Wigner functions. E.g., the latter can 
be given by |j| 



WM(p)=N<%- 2 W° 



JA-l(V=8i^) 



(0*)(A-l}/2 



(40) 



where Jx-i(x) is the Bessel function and is the nor- 
malization constant. 

To find a physical insight of the steady-state density 
matrix (|38|l . we rewrite it as the incoherent superposition 
(PJ): 



: Po|V'o)(V'o| +Pl\lpl)(i>l\ +P2\lp2)(lp2 



(41) 



of the orthogonal pure states \ipk) with the weight factors 
Pk, which can be found numerically. Applying the eigen- 
value method for the steady state we have obtained all 
the spectrum of the final state. Wigner functions of the 
first three most important components are displayed in 
figure 11. The weight coefficients of the strongest com- 
ponents are in the sequence of po as 0.928, p\ as 0.068 
and pi w 0.004. The first component |i/>o)j apart from 
the second and third components, is strongly nonclassi- 
cal with its quadrature noise of about So ~ 0.60, while 
the compound steady state has the quadrature noise level 
just about S pa 0.72. To verify the conclusion that the 
density matrices p(t) for different initial conditions evolve 
into the same steady state described by the density ma- 
trix H38|) . we calculate the Bures metric and the quan- 
tum relative entropy (the generalized Kullback-Leibler 
distance) between p(t) and p ss defined as follows ^(| 



D B {p(t)\\p BS } 
D KL {p(t)\\p ss } 



2Tr{[, 



2 },(42) 



Tr{p{t)[lnp{t)-lnp ss ]}, (43) 



respectively. The square of trace in l|42|) is the so-called 
Uhlmann transition probability or the Jozsa fidelity for 
mixed states Figure 10 clearly shows that the Bures 
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metric and quantum relative entropy approach zero for 
the three different initial fields including coherent states, 
their superpositions and the Fock states. 

In the next section we will show that the steady state 
can be well described as a nonclassical Gaussian state at 
least for high-pump intensities. 



V. LINEARIZED APPROXIMATION OF THE 
STEADY-STATE SOLUTION 

Now we will give a linearized approximation of the 
steady-state solution of the Heisenberg equation © . We 
start with the classical approximation for strong fields, 
when the quantum noise can be neglected and the signal 
mode can be described by a complex amplitude. Heisen- 
berg equation |JS} gives the following ordinary differential 
equation 



do; n 

— = p - 2iG\a\ 2 a - j a 



(44) 



for the complex amplitude a. Minor differences between 
the exact quantum mean amplitude (a (t)) and the clas- 
sical amplitude a (t) for the coherent initial state \a = 3) 
and for the resonator parameters p = 5, G = 0.2 and 
7o = 1 can be seen in the phase-space trajectory depicted 
in figure 9. 

If we are interested in the steady-state solution only, 
we may replace the left hand side of equation (|44|l by zero 
and solve it. The complex equation can be transformed 
into the real equation 



(7o 



■4G 2 H 4 ) la! 2 , 



which is a cubic equation in intensity 



(45) 
A simple 

analysis shows that the amplitude |a| of the steady-state 
solution is a monotonic function of the pump parame- 
ter \p\ and so equation (|45|l has a unique solution and 
no threshold can be observed. Since the explicit form of 
the solution is too complex to find some relevant physics 
in it, we will not present it here. We note only that 
for the resonator parameters p — 5, G — 0.2 and 70 = 1 
assumed in the examples investigated in the previous sec- 
tion one can find the exact classical steady-state solution 
\a\ 2 — 5 and a = 1 — 2i. In particular, (i) for p = 0, 
equation (|45|1 has a trivial solution of a = 0, (ii) for 
weak pump \p\ 2 <C 7q/2G, the amplitude grows linearly 
with the pump intensity as a ~ p/70 and (iii) for very 
strong pump \p\ 2 3> 7q/2G, the amplitude grows with 
the third root of the pump intensity 



P_ 

Vi 



2G 



1/3 



(46) 



To describe the noise properties of the signal mode we 
will search a solution to equation as a sum of two 
components a = a + A, where a is the strong classi- 
cal complex amplitude and A is the weak quantum noise 



operator. By inserting it into equation © and after lin- 
earization one gets the following equation for the noise 
operator 



(47) 



where the complex coefficients 7 and 5 depend on time 
through the amplitude a and are given by the relations 



7 = 7o + 4iG|a| :! 



2\Ga 2 



(48) 



The solution of equation l|47|l is a Gaussian state, defined 
by ljA2fl in Appendix A, exhibiting purely nonclassical 
properties, thus referred to as the NCGS. The time evo- 
lution of the NCGS parameters B and G is described by 
the two differential equations 



B = -2( 7 + 7*)B- (S*C 
C = -S(l + 2B)-2jC. 



SC*) 



(49) 



From them one can obtain numerically the linearized so- 
lution of the Heisenberg equation @, but it is mean- 
ingless when we know the exact quantum solution. In- 
stead of that we will focus on the steady-state solution 
only. Since I7I > \S\, a stationary solution of the operator 
equation (|47H exists. It is the NCGS with the following 
parameters |l8j 



2B = 



\5f 



\l\ 2 -\5f 



2C 



S-y* 



h\ 2 -\6\< 



(50) 



Since we know all the three parameters a, B and G of the 
NCGS, we can simply estimate the other parameters of 
the steady state, where the corresponding formulas are 
given in Appendix A. For example, for the quadrature 
noise squeezing parameter (|A7|) one can derive a simple 
formula 



S = 



\l\ 



h\ + \s\ 



< 1 



(51) 




demonstrating that the steady state is squeezed. The 
parameter x, defined by (|A4|) . is explicitly given by 



(52) 



Numerically, for the resonator parameters p = 5, G = 0.2 
and 70 = 1, we get a = 1 — 2i, 7 = 1 + 4i, S = 1. 6 — 1. 
2i, S w 0.673, F « 0.711, x w 0.072, p « 0.933, p 1 w 
0.062 and P2 ~ 0.004. From formula (|A6|) we estimate 
the linear entropy of the steady state as L w 0.126 and 
from (|A5|) we estimate the von Neumann entropy as E rj 
0.263. Surprisingly, these estimations match the exact 
values already for relatively low pump intensities. For 
higher pump intensities the agreement is even better. 

Note that in the strong pump approximation, a very 
simple relation holds I7I = 2\8\, which gives crude es- 
timations of S — 2/3 and F — 2/3. Since x = 



9 




cold reservoir. 
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FIG. 12: Exact quantum weight factors Pk of the steady state 
(thick curve with squares) vs the Gaussian state approxima- 
tion given by <A3t (thin curve with circles) 



|v3 — \ w 0.0774, the linear entropy is estimated as 
L = l-y/3/2 « 0.134 and the entropy as E « 0.278. We 
can also estimate the weight coefficients of steady state 
decomposition as po » 0.928, p\ w 0.067, pi w 0.005. So, 
even the crude estimations well match the exact values 
obtained numerically in the previous section. 

In figure 12 we have compared the exact quantum 
weight factors pk of the steady state and the approximate 
Gaussian weight factors given by l|A3|l with x w 0.072 cal- 
culated from (|A4|) for the resonator parameters G = 0.2, 
7o = 1 and p — 5. This graphical representation shows 
the validity of our approximation. 



VI. CONCLUSIONS 

We studied the evolution of the quantum noise and 
mixedness of a dissipative non-linear oscillator, described 
by a Kerr non-linear oscillator, pumped by a classical ex- 
ternal field. We quantified the quantum noise by the 
quadrature squeezing parameter and the Fano factor, 
while the mixedness by the quantum von Neumann and 
linear entropies. Dissipation was described in the stan- 
dard Heisenberg-Langevin approach by coupling the sys- 
tem to zero-temperature reservoir. We demonstrated 
that initial pure states, including the Fock states, co- 
herent states, and their finite superpositions, exhibited 
fast decoherence and evolved into the same steady mixed 
state being well approximated by a nonclassical Gaus- 
sian state as verified by the Bures metric and the quan- 
tum relative entropy. We presented analytical formulas 
and numerical results for the steady state parameters to 
show that the state exhibits quadrature noise squeezing 
and sub-Poissonian statistics if the thermal reservoir is 
cold. We observed, especially for the initial nonclassical 
states including the Fock states or the finite superposi- 
tions of coherent states, a rapid increase in the mixedness 
at the very beginning of the evolution, and then for longer 
times a fall in the mixedness to the same low value of the 
asymptotic nonclassical Gaussian state in the case of the 



APPENDIX A: GAUSSIAN STATES 

The Gaussian state is defined as a state with Gaussian 
quasidistribution and Gaussian characteristic functions 
(see, e.g., 0j^3H3)i which can be identified completely 
by the three parameters 

a={a), B = (tfa)-(tf)(a), C = (a 2 )-(a) 2 . (Al) 

Thus the Gaussian state can be defined by its Cahill- 
Glauber s-parametrizcd quasidistribution function 



W (fl) (/3) = 



1 



exp 



1 - s + 2B 
2Kl 



1/3- «| S 



— (f3 - af + c.c. 



(A2) 



where K s = (1/2 - s/2 + B) 2 — \C\ 2 . The Gaussian state 
is a natural generalization of coherent [B = C = 0) and 
chaotic (C = 0) states. Note that in a special case of 
s = 1, the Glauber-Sudarshan function P(0) = W^\0) 
exists if Ki > 0, otherwise the state described by (|A2|> 
is nonclassical, and thus referred to as the NCGS. The 
Gaussian states are mixed states with the weight coeffi- 
cients 



Pk 



(1 



,i+fe ' 



(A3) 



where 




(A4) 



The von Neumann entropy of the Gaussian state 
reads as fl3l 



In- 



(A5) 



(i + x) 1 

Similarly, it is easy to show that the purity parameter is 
P=(l + 2xy\ (A6) 
while the quadrature noise squeezing parameter reads as 
S=l + 2(B-\C\) (A7) 



10 



and its Fano factor is where the parameters a, B and C are defined by i|Al|l . 

Note that the squeezing occurs if B < \C\ or equivalcntly 

f V"V* —I— C~** C\i 

F = 1 + 2B + (A8) i 1 < 0' so an y NCGS exhibits quadrature squeezing. 

|a| 2 +i? 
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